第 2 章:和式

记号

求和的符号有两种形式

第一种是确定界限的形式,也叫封闭形式,例如:k=1nak

第二种叫做一般形式,就是把一个或者多个条件写在 符号的下面,例如刚刚的例子可以写成 1knak

和式和递归式的转化

和式和递归式之间可以相互转化,

和式转化成递归式

Sn=an 可以转化为 Sn=Sn1+an,S1=a1,后面这一项就是一个递归式子

例如:我们这里有一个递归式

R0=αRn=Rn1+β+γn,n>0

这个可以写成通解的形式:Rn=A(n)α+B(n)β+C(n)γ

由成套方法得出 A(n)=1,B(n)=n,C(n)=(n2+n)/2

我们需要计算 k=0n(a+bk),只需要把 a=α,a=β,b=γ 代入即可

递归式转化成和式

T0=0Tn=2Tn1+1,n>0

两边除以 2n 得到

T0/20=0Tn/2n=Tn1/2n1+1/2n,n>0

Sn=Tn/2n

S0=0Sn=Sn1+1/2n,n>0

由此得到 Sn=k=1n12k=1(12)n,所以 Tn=2n1

和式的变换法则

经典的三条性质,注意这里的 K 是一个集合

kKcak=ckKak(分配律)kK(ak+bk)=kKak+kKbk(结合律)kKak=p(k)Kap(k)(交换律)

注意 p(k) 是一个双射函数,也可以理解为一个排列

例如:K={1,2,3,4} 那么,p(k)={4,2,3,1} 就是一个排列

当然 p(k) 可以多出几个元素例如上面那个例子 p(k)={4,2,3,1,5,6} 也是可以的,5,6 不在 K

接下来看一个例子:S=0kn(a+bk)

利用交换率,用 nk 代替 k,得到

S=0nkn(a+b(nk))=0kn(a+bnbk)

利用结合律,把这两个方程相加,得到

2S=0kn(a+bk+a+bnbk)=0kn(2a+bn)=(n+1)(2a+bn)

所以

S=(n+1)(2a+bn)2

扰动法

和式中有一种非常厉害的方法叫扰动法,就是把一项从和式中去除出去,然后尝试把剩下的项变成 Sn 的形式,从而解出 Sn

例如:Sn=0knaxk,使用扰动法的基本套路

Sn+axn+1=ax0+1kn+1axk=ax0+0knaxk+1=ax0+xSn

等式两边同时出现了 Sn,解出 Sn=aaxn+11x

考虑另外一个例子:Sn=0knk2k,考虑扰动法

Sn+(n+1)2n+1=0kn(k+1)2k+1=20knk2k+0kn2k+1=2Sn+2n+22

解出 Sn=2n+1(n1)+2

设未知数 x,能求出 Sn=0knkxk 的通解:

Sn=xxn+1(n+1)+nxn+2(x1)2,x1

我们也可以利用求导的方法,求出 Sn=0knkxk 的通解:

已知:k=0nxk=1xn+11x,两边同时求导

k=0nkxk1=(1x)((n+1)xn)+1xn+1(1x)2=1(n+1)xn+nxn+1(1x)2

两边同时乘以 x,得到

k=0nkxk=xxn+1(n+1)+nxn+2(x1)2

就能得到和上面扰动法一样的结果了

多重和式

多重和式对应于连续函数的积分,可以利用积分的一些思维来思考多重和式

对于自变量相互无关的情况很好理解

1j,k3ajbk=a1b1+a1b2+a1b3+a2b1+a2b2+a2b3+a3b1+a3b2+a3b3=a1(b1+b2+b3)+a2(b1+b2+b3)+a3(b1+b2+b3)=(a1+a2+a3)(b1+b2+b3)=(1j3aj)(1k3bk)

一个多重和式可以转化成两个和式相乘的形式,前提是 J=K,其中 jJ,kK,这一点在级数相乘的时候极为有效

在具体求多重和式的时候,往往是一层一层往外求,理论上来说,先固定住哪个自变量都是无所谓的,所以我们可以交换求和符号的位置

P(j,k) 返回的是 j,k 是否满足某种性质,那么有:

jJkKaj,k[P(j,k)]=kKjJaj,k[P(j,k)]

讨论自变量之间有某些限制

一种比较常见的方式:1jkn,我们可以通过画图找出需要求和的元素

image.png|512

可以看到对角线以及右上角的那一块都是需要求和的

另外一种常见的方式 1j,kn,j+k=n,画出图来是按对角线求和

image.png

这种求和方法有一种自己的名字,叫卷积(Convolution)

例一:求 1jknajak

image.png

可以看到要求部分是右上三角形

由于图标的对称性,可以得到 S左下=S右上

由容斥原理:S左下+S右上=S全部S对角线

于是就能得我们想要的答案了

S右上=1jknajak=12((k=1nak)2+k=1nak2)

例二:求 S=1j<kn(akaj)(bkbj)

交换 j,k 仍然由对称性:

S=1k<jn(ajak)(bjbk)=1k<jn(akaj)(bkbj)

所以得到答案

2S=1j,kn(ajak)(bjbk)1j=kn(ajak)(bjbk)

显然,第二个和式等于 0,把第一个和式展开

2S=1j,knajbj1j,knajbk1j,knakbj+1j,knakbk=2n1knakbk2(k=1nak)(k=1nbk)

所以就得出了 S=n1knakbk(k=1nak)(k=1nbk)

这两个不等式叫做 切比雪夫单调不等式

例三

Sn=1j<kn1kj

我们可以用 kj 替换 j,然后固定 k

Sn=1kn1jk1kj=1kn0jk11j=1knHk1=0kn1Hk

其中 Hk 是调和级数,我们对调和级数有公式,但是对于调和级数求和没有公式,说明固定 k 走不通考虑固定 j

Sn=1j<n j<kn1kj=1jn j<k+jn1k+jj=1jn 0<knj1k=1jnHnj=0jn1Hj

还是回到了调和级数求和的问题,考虑按对角线求和

Sn=1j<n j<kn1kj=1jn j<k+jn1k+jj=1jn 0<knj1k=0<k<n 1jnk1k=0<k<nnkk=0<k<n(nk1)=n1kn1knn(n1)=nHnn

我们还可以从几何的层面来理解这个式子:

image.png|421

刚开始的几次我们按照行和列求和,都会得到一个调和级数求和的式子,最后一种方式是按照对角线,这里得到 31+22+13

一般性方法

这一节基本上就是总结性质的,把前面介绍的方法综合运用一下

要求一个立方和 Sn=k=0nk2, n0

方法0:查找公式

OEIS 能帮我们很好的查找公式

得到的答案是 Sn=n(n+1)(2n+1)6

方法1:数学归纳法

已知:Sn=n(n+12)(n+1)3

显然 S0=0 成立,假设 n>0 时,Sn1=(n1)(n12)n3 成立,有:Sn=Sn1+S

3Sn=(n1)(n12)(n)+3n2=n3+32n2+12n=n(n+1)(n+12)

所以 Sn=n(n+12)(n+1)3 对所有 n0 都成立

方法2:扰动法

Sn+(n+1)2=k=0n(k+1)2=k=0n(k2+2k+1)=k=0nk2+2k=0nk+k=0n1=Sn+2k=0nk+n+1

惊人的发现,我们需要求的 Sn 被消掉了,留下一个 k=0nk=(n+1)2(n+1)2=n(n+1)2

我们猜想能不能用立方和来把平方和求出来,设 Tn=k=0nk3

Tn+(n+1)3=k=0n(k+1)3=k=0n(k3+3k2+3k+1)=k=0nk3+3k=0nk2+3k=0nk+k=0n1=Tn+3Sn+3n(n+1)2+n+1

Tn 被消掉了,留下我们需要求的 Sn,这样 3Sn=(n+1)33n(n+1)2(n+1)

所以 Sn=n(n+1)(n+12)3

方法3:建立成套方法

还是老规矩,建立成套方法,设:

R0=α,n=0Rn=Rn1+β+γn+δn2,n>0

解的一般形式就是:Rn=A(n)α+B(n)β+C(n)γ+D(n)δ

我们已知,δ=0 时,A(n)=1,B(n)=n,C(n)=(n2+n)/2

现在要求 D(n) 和他们的关系,设 Rn=n3,有 n3=(n1)3+β+γn+δn2

可得:α=0β=1γ=3δ=3,代入:3D(n)3C(n)+B(n)=n3

就把 D(n) 接出来了,D(n)=n33n(n+1)2+n=n(n+12)(n+1)3

我们需要求的和式 Sn 转化成递归式子,只需要令 α=β=γ=0 以及 δ=1

Sn=D(n)

方法4:用积分替换和式

我们已知 0nx2dx=n33,显然,和式和积分之间的差距并不差,只是差一个误差

我们设这个误差为 En=Snn33,通过 Sn 的递归式,我们能得出 Sn=Sn1+n2

En=Snn33=Sn1+n2n33=En1+(n1)33+n2n33=En1+n13

于是,我们就得到了 En 的递推式,En=En1+n13E0=0

很容易算出 En 的通项公式为 En=n(n+1)2n3,所以 Sn=En+n33=n(n+1)(n+12)3

方法5:展开和放缩

这种方法极具技巧性,把一个一重和式转化成二重和式

Sn=k=0nk2=1jknk=1jn jknk=1jn(j+n2)(nj+1)=121jn(n(n+1)+jj2)=12n2(n+1)+14n(n+1)12Sn

这样等式两边都出现了 Sn,解出 Sn=n(n+1)(2n+1)6

有限微积分

我们类似于积分和微分的思维,在离散域上定义有限微积分,定义差分算子 Δf(x)(f(x+h)f(x))/hh=1 时的值

还需要定义两种特殊的次幂:

我们为什么要定义这样奇特的形式呢,因为这样求差分的时候会和微分有很多共同点

Δxm=(x+1)mxm=(x+1)x(x1)(xm+2)x(x1)(x2)(xm+1)=(x+1(xm+1))x(x1)(x2)(xm+2)=mx(x1)(x2)(xm+2)=mxm1

有限微积分也存在类似于微积分基本定理的东西

g(x)=Δf(x)g(x)δx=f(x)+C

这里的 g(x)δxg(x) 的不定和式,这里的 C 可以是周期为 1 的任意函数

大名鼎鼎的牛顿-莱布尼茨公式在有限微积分中也有体现 abg(x)δx=f(x)|ab=f(b)f(a)

思考 abg(x)δx=f(x)|ab=f(b)f(a) 当,上下限相同时:

aag(x)δx=f(x)|aa=f(a)f(a)=0

当上下限差 1 时:

aa+1g(x)δx=f(x)|aa+1=f(a+1)f(a)=Δf(x)=g(x)

着预示着和传统的求和符号有点不同,可以理解为 abg(x)δx=k=ab1g(k)=ak<bg(k),ba

也就是说,如果存在能找到一个 f(x) 使得 f(x+1)f(x)=g(x),那么理论上来说就可以求 g(x) 这类的和了

我们尝试去求类似于原函数的东西 f(x)

比如:g(x)=xmf(x)=xm+1m+1

我们用前面的例子,求 0k<nk2 已知 k2=k2+k1,所以

0k<nk2=n33+n22=n(n1)(n2+32)=13n(n12)(n1)

利用 n+1 替换 n 就可以得到 Sn

尝试继续推广无限微积分,考虑下降次幂为负数的情况,xm, m<0

观察:

于是找规律得出 x1=1x+1,x2=1(x+1)(x+2)xm=1(x+1)(x+2)(x+m)

通常的幂法则:xm+n=xmxn,推广到有限微积分就变成了 xm+n=xm(xm)n

现在确认下降幂的差分性质在 m<0 是否成立也就是 Δxm=mxm1

如果 m=2,有

Δx2=(x+1)2x2=1(x+2)(x+3)1(x+1)(x+2)=(x+1)(x+3)(x+1)(x+2)(x+3)=2(x+1)(x+2)(x+3)=2x3

说明求和性质对 m<0 也成立,即 abxmδx=xm+1m+1|ab, m1

现在考虑 m=1 的时候,对于微积分,我们又 abx1dx=lnx|ab

我们需要在有限微积分中找一个类似于 lnx 的函数 ,这个函数的差分为 1x+1

很显然能得到这个函数就是 Hx,于是就得到了求和的完整形式

abxmδx={xm+1m+1|ab,m1,Hx|ab,m=1.

类似的,我们需要找一个 ex 类似物,根据定义 d(ex)=ex,所以有 f(x+1)f(x)=f(x)f(x+1)=2f(x),即 f(x)=2x

cx 的差分也相当简单,即对任意的 c

Δ(cx)=cx+1cx=(c1)cx.

那么,c1 时,cx 的原函数就是 cx/(c1)

两个相乘的函数求微分:d(uv)=udv+vdu,两边同时积分得到分部积分法的公式:udv=uvvdu

有限微积分也类似

Δ(u(x)v(x))=u(x+1)v(x+1)u(x)v(x)=u(x+1)v(x+1)u(x)v(x+1)+u(x)v(x+1)u(x)v(x)=u(x)Δv(x)+v(x+1)Δu(x).

这里的 v(x+1) 看着很烦,所以定义一个位移算子 Ef(x)=f(x+1)

所以乘积差分法则可以写为:

Δ(uv)=uΔv+EvΔu

两边求和得:

uΔv=uvEvΔu

有限微积分和无限微积分对应着离散和连续,有很多共通点,在下列表中给出

image.png|789

来看几个例子:

例一:

k=0nk2k=0n+1x2xδx=x2x2x+1|0n+1=((n+1)2n+12n+2)(0×2021)=(n1)2n+1+2

例二:0k<nkHk

先求原函数,根据分部积分法则:

xHxδx=x22Hx(x+1)2xx1δx=x22Hx12x1δx=x22Hxx24+C

然后把上下限代入:

0k<nkHk=0nxHxδx=n22(Hn12)

无限和式

引用大乌拉的一句话:

我们对无穷的东西几乎一无所知

这一章书上举了几个具体的例子,都在表达无穷和式似乎没有我们想象的这么简单